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A detailed analysis of the November f5, 2006 data release (Bradac et al., 2006; Clowe et al., 
2006a, b,c) X-ray surface density S-map and the strong and weak gravitational lensing conver- 
gence /t-map for the Bullet Cluster fE0657-558 is performed and the results are compared with 
the predictions of a modified gravity (MOG) and dark matter. Our surface density S-model is 
computed using a King /3-model density, and a mass profile of the main cluster and an isothermal 
temperature profile are determined by the MOG. We find that the main cluster thermal profile 
is nearly isothermal. The MOG prediction of the isothermal temperature of the main cluster is 
T = 15.5 ± 3.9 keV, in good agreement with the experimental value T = 14.8^^ keV. Excellent 
fits to the two-dimensional convergence K-map data are obtained without non-baryonic dark mat- 
ter, accounting for the 8<r spatial offset between the S-map and the K-map reported in Clowe et al. 
(2006a). The MOG prediction for the K-map results in two baryonic components distributed across 
the Bullet Cluster fE0657-558 with averaged mass-fraction of 83% intracluster medium (ICM) 
gas and f 7% galaxies. Conversely, the Newtonian dark matter K-model has on average 76% dark 
matter (neglecting the indeterminant contribution due to the galaxies) and 24% ICM gas for a 
baryon to dark matter mass-fraction of 0.32, a statistically significant result when compared to 
the predicted A-CDM cosmological baryon mass- fraction of 0.176±o' o il (Spergel et al, 2007). 
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I. INTRODUCTION 

A. The Question of Missing Mass 

Galaxy cluster masses have been known to require some form of energy density that makes its presence felt only by 
its gravitational effects since Zwicky (1933) analysed the velocity dispersion for the Coma cluster. Application of the 
Newtonian 1/r 2 gravitational force law inevitably points to the question of missing mass, and may be explained by 
dark matter (Oort, 1932). The amount of non-baryonic dark matter required to maintain consistency with Newtonian 
physics increases as the mass scale increases so that the mass to light ratio of clusters of galaxies exceeds the mass to 
light ratio for individual galaxies by as much as a factor of ~ 6, which exceeds the mass to light ratio for the luminous 
inner core of galaxies by as much as a factor of ~ 10. In clusters of galaxies, the dark matter paradigm leads to a mass 
to light ratio as much as 300 Mq/Lq. In this scenario, non-baryonic dark matter dominates over baryons outside the 
cores of galaxies (by » 80 — 90%). 

Dark matter has dominated cosmology for the last five decades, although the search for dark matter has to 
this point come up empty. Regardless, one of the triumphs of cosmology has been the precise determination of 
the standard (power-law flat, A-CDM) cosmological model parameters. The highly anticipated third year results 
from the WMAP team have determined the cosmological baryon mass-fraction (to non-baryonic dark matter) to be 
0.176^q'q^2 (Spergel et at, 2007). This ratio may be inverted, so that for every gram of baryonic matter, there are 
5.68 grams of non-baryonic dark matter - at least on cosmological scales. There seems to be no evidence of dark 
matter on the scale of the solar system, and the cores of galaxies also seem to be devoid of dark matter. 

Galaxy clusters and superclusters are the largest virialized (gravitationally bound) objects in the Universe and 
make ideal laboratories for gravitational physicists. The data come from three sources: 

i. X-ray imaging of the hot intracluster medium (ICM), 

ii. Hubble, Spitzer and Magellan telescope images of the galaxies comprising the clusters, 

iii. strong and weak gravitational lensing surveys which may be used to calculate the mass distribution projected 
onto the sky (within a particular theory of gravity). 

The alternative to the dark matter paradigm is to modify the Newtonian 1/r 2 gravitational force law so that the 
ordinary (visible) baryonic matter accounts for the observed gravitational effect. An analysis of the Bullet Cluster 
1E0657-558 surface density S-map and convergence K-map data by Angus et al. (2006, 2007) based on Milgroms Modi- 
fied Newtonian Dynamics (MOND) model (Milgrom, 1983; Sanders and McGaugh, 2002) and Bekenstein's relativistic 
version of MOND (Bekcnstcin, 2004) failed to fit the data without dark matter. More recently, further evidence that 
MOND needs dark matter in weak lensing of clusters has been obtained by Takahashi and Chiba (2007). Problems 
with fitting X-ray temperature profiles with Milgroms MOND model without dark matter were shown in Aguirre et al. 
(2001); Brownstein and Moffat (2006a); Pointecouteau and Silk (2005); Sanders (2006). Neutrino matter with an elec- 
tron neutrino mass m„ ~ 2 eV can fit the bullet cluster data (Angus et al, 2006, 2007; Sanders, 2006). This mass is 
at the upper bound obtained from observations. The Karlsruhe Tritium Neutrino (KATRIN) experiment will be able 
to falsify 2 eV electron neutrinos at 95% confidence level within months of taking data in 2009. 

The theory of modified gravity - or MOG model - based upon a covariant generalization of Einstein's theory with 
auxiliary (gravitational) fields in addition to the metric was proposed in Moffat (2005, 2006b) including metric skew- 
tensor gravity (MSTG) theory and scalar-tensor-vector gravity (STVG) theory. Both versions of MOG, MSTG and 
STVG, modify the Newtonian 1/r 2 gravitational force law in the same way so that it is valid at small distances, say 
at terrestrial scales. 

Brownstein and Moffat (2006b) applied MOG to the question of galaxy rotation curves, and presented the fits to a 
large sample of over 100 low surface brightness (LSB), high surface brightness (HSB) and dwarf galaxies. Each galaxy 
rotation curve was fit without dark matter using only the available photometric data (stellar matter and visible gas) 
and alternatively a two-parameter mass distribution model which made no assumption regarding the mass to light 
ratio. The results were compared to MOND and were nearly indistinguishably right out to the edge of the rotation 
curve data, where MOND predicts a forever flat rotation curve, but MOG predicts an eventual return to the familiar 
1/r 2 gravitational force law. The mass to light ratio varied between 2 — 5 Mq/Lq across the sample of 101 galaxies 
in contradiction to the dark matter paradigm which predicts a mass to light ratio typically as high as 50 Mq/Lq. 

In a sequel, Brownstein and Moffat (2006a) applied MOG to the question of galaxy cluster masses, and presented 
the fits to a large sample of over 100 X-ray galaxy clusters of temperatures ranging from 0.52 keV (6 million kelvin) 
to 13.29 keV (150 million kelvin). For each of the 106 galaxy clusters, the MOG provided a parameter-free prediction 
for the ICM gas mass profile, which reasonably matched the X-ray observations (King /3-model) for the same sample 
compiled by Reiprich (2001) and Reiprich and Bohringcr (2002). The MOND predictions were presented for each 
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galaxy cluster, but failed to fit the data. The Newtonian dark matter result outweighed the visible ICM gas mass 
profiles by an order of magnitude. 

In the solar system, the Doppler data from the Pioneer 10 and 11 spacecraft suggest a deviation from the Newtonian 
1/r 2 gravitational force law beyond Saturn's orbit. Brownstein and Moffat (2006c) applied MOG to fit the available 
anomalous acceleration data (Nieto and Anderson, 2005) for the Pioneer 10/11 spacecraft. The solution showed 
a remarkably low variance of residuals corresponding to a reduced x 2 P er degree of freedom of 0.42 signalling a 
good fit. The magnitude of the satellite acceleration exceeds the MOND critical acceleration, negating the MOND 
solution (Sanders, 2006). The dark matter paradigm is severely limited within the solar system by stability issues 
of the sun, and precision gravitational experiments including satellite, lunar laser ranging, and measurements of the 
Gaussian gravitational constant and Kepler's law of planetary motion. Without an actual theory of dark matter, no 
attempt to fit the Pioneer anomaly with dark matter has been suggested. Remarkably, MOG provides a closely fit 
solution to the Pioneer 10/11 anomaly and is consistent with the accurate equivalence principle, all current satellite, 
laser ranging observations for the inner planets, and the precession of perihelion for all of the planets. 

A fit to the acoustical wave peaks observed in the cosmic microwave background (CMB) data using MOG has been 
achieved without dark matter. Moreover, a possible explanation for the accelerated expansion of the Universe has 
been obtained in MOG (Moffat, 2007). 

Presently, on both an empirical and theoretical level, MOG is the most successful alternative to dark matter. The 
successful application of MOG across scales ranging from clusters of galaxies (Megaparsecs) to HSB, LSB and dwarf 
galaxies (kiloparsecs) , to the solar system (AU's) provides a clue to the question of missing mass. The apparent 
necessity of the dark matter paradigm may be an artifact of applying the Newtonian 1/r 2 gravitational force law to 
scales where it is not valid, where a theory such as MOG takes over. The "excess gravity" that MOG accounts for 
may have nothing to do with the hypothesized missing mass of dark matter. But how can we distinguish the two? 
In most observable systems, gravity creates a central potential, where the baryon density is naturally the highest. 
So in most situations, the matter which is creating the gravity potential occupies the same volume as the visible 
matter. Clowe et al. (2006c) describes this as a degeneracy between whether gravity comes from dark matter, or 
from the observed baryonic mass of the hot ICM and visible galaxies where the excess gravity is due to MOG. This 
degeneracy may be split by examining a system that is out of steady state, where there is spatial separation between 
the hot ICM and visible galaxies. This is precisely the case in galaxy cluster mergers: the galaxies will experience 
a different gravitational potential created by the hot ICM than if they were concentrated at the center of the ICM. 
Moffat (2006a) considered the possibility that MOG may provide the explanation of the recently reported "extra 
gravity" without non-baryonic dark matter which has so far been interpreted as direct evidence of dark matter. The 
research presented here addresses the full-sky data product for the Bullet Cluster 1E0657-558, recently released to 
the public (Clowe et al, 2006b). 




FIG. 1: 

False colour image of Bullet Cluster 1E0657-558. 

The surface density E-map reconstructed from 
X-ray imaging observations is shown in red 
and the convergence K-map as reconstructed 
from strong and weak gravitational lensing ob- 
servations is shown in blue. Image pro- 
vided courtesy of Chandra X-ray Observatory. 
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B. The Latest Results from the Bullet Cluster 1E0657-558 

• The Chandra Peer Review has declared the Bullet Cluster 1E0657-558 to be the most interesting cluster in 
the sky. This system, located at a redshift z = 0.296 has the highest X-ray luminosity and temperature 
(T = 14.1 ± 0.2 keV ~ 1.65 x 10 s K), and demonstrates a spectacular merger in the plane of the sky exhibiting a 
supersonic shock front, with Mach number as high as 3.0 ± 0.4 (Markevitch, 2006). The Bullet Cluster 1E0657- 
558 has provided a rich dataset in the X-ray spectrum which has been modeled to high precision. From the 
extra-long 5.2 x 10 5 s Chandra space satellite X-ray image, the surface mass density, y), was reconstructed 
providing a high resolution map of the ICM gas (Clowe et ai, 2006c). The S-map, shown in a false colour 
composite map (in red) in Figure 1 is the result of a normalized geometric mass model based upon a 16' x 16' 
field in the plane of the sky that covers the entire cluster and is composed of a square grid of 185 x 185 pixels 
(~ 8000 data-points) 1 . Our analysis of the E-map provides first published results for the King /3-model density, 
as shown in Table 2, and mass profiles of the main cluster and the isothermal temperature profile as determined 
by MOG, as shown in Figure 9. 

• Based on observations made with the NASA/ESA Hubble Space Telescope, the Spitzer Space Telescope and 
with the 6.5 meter Magellan Telescopes, Bradac et al. (2006); Clowe et al. (2006a,c) report on a combined strong 
and weak gravitational lensing survey used to reconstruct a high-resolution, absolutely calibrated convergence 
re-map of the region of sky surrounding Bullet Cluster 1E0657-558, without assumptions on the underlying 
gravitational potential. The re-map is shown in the false colour composite map (in blue) in Figure 1. The 
gravitational lensing reconstruction of the convergence map is a remarkable result, considering it is based on a 
catalog of strong and weak lensing events and relies upon a thorough understanding of the distances involved - 
ranging from the redshift of the Bullet Cluster 1E0657-558 (z = 0.296) which puts it at a distance of the order 
of one million parsecs away. Additionally, the typical angular diameter distances to the lensing event sources 
(z ~ 0.8 to z ~ 1.0) are several million parsecs distant! This is perhaps the greatest source of error in the re-map 
which limits its precision. Regardless, we are able to learn much about the convergence map and its peaks. 

Both the E-map and the re-map are two-dimensional distributions based upon line-of-sight integration. We are 
fortunate, indeed, that the Bullet Cluster 1E0657-558 is not only one of the hottest, most supersonic, most massive 
cluster mergers seen, but the plane of the merger is aligned with our sky! As exhibited in Figure 1, the latest results 
from the Bullet Cluster 1E0657-558 show, beyond a shadow of doubt, that the S-map, which is a direct measure of 
the hot ICM gas, is offset from the re-map, which is a direct measure of the curvature (convergence) of space-time. 
The fact that the re-map is centered on the galaxies, and not on the ICM gas mass is certainly either evidence of 
"missing mass", as in the case of the dark matter paradigm, or "extra gravity", as in the case of MOG. Clowe et al. 
(2006c) states 

One would expect that this (the offset S- and re-peaks) indicates that dark matter must be present 
regardless of the gravitational force law, but in some alternative gravity models, the multiple peaks can 
alter the lensing surface potential so that the strength of the peaks is no longer directly related to the 
matter density in them. As such, all of the alternative gravity models have to be tested individually 
against the observations. 

The data from the Bullet Cluster 1E0657-558 provides a laboratory of the greatest scale, where the degeneracy between 
"missing mass" and "extra gravity" may be distinguished. We demonstrate that in MOG, the convergence re-map 
correctly accounts for all of the baryons in each of the main and subclusters, including all of the galaxies in the regions 
near the main central dominant (cD) and the subcluster's brightest central galaxy (BCG), without non-baryonic dark 
matter. 



C. How Modified Gravity may account for the Bullet Cluster 1E0657-558 evidence 

We will show how MOG may account for the Bullet Cluster 1E0657-558 evidence, without dominant dark matter, by 
deriving the modifications to the gravitational lensing equations from MOG. Concurrently, we will provide comparisons 
to the equivalent Einstein-Newton results utilizing dark matter to explain the missing mass. The paper is divided as 
follows: 



1 Technical details in Markevitch, M. et al. 2007, in preparation. 
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Section II is dedicated to the theory used to perform all of the derivations and numerical computations and is 
separated into three pieces: Section II. A presents the Poisson equations in MOG for a non-spherical distribution 
of matter and the corresponding derivation of the acceleration law and the dimensionless gravitational coupling, Q. 
Section II. B presents the King /3-model density profile, p. Section II. C presents the derivation of the weighted surface 
mass density, E = J Qp, from the convergence re — E/E c . The effect of the dimensionless gravitational coupling, Q, is 
to carry more weight away from the center of the system. If the galaxies occur away from the center of the ICM gas, 
as in the Bullet Cluster 1E0657-558, their contribution to the re-map will be weighted by as much as a factor of 6 as 
shown in Figure 7. 

Section III is dedicated to the surface density E-map from the X-ray imaging observations of the Bullet Cluster 
1E0657-558 from the November 15, 2006 data release E-map (Clowe et at, 2006b). Section III. A presents a visualiza- 
tion of the E-map and our low x 2 best- fit King /3-model (neglecting the subcluster). Section III.B presents a determi- 
nation of Q for the Bullet Cluster 1E0657-558 based upon the (> 100) galaxy cluster survey of Brownstein and Moffat 
(2006a). Section III.C presents the cylindrical mass profile, J E, about the main cluster E-map peak. Section III.D 
presents the isothermal spherical mass profile from which we have derived a parameter-free (unique) prediction for 
the X-ray temperature of the Bullet Cluster 1E0657-558 of T = 15.5 ± 3.9 keV which agrees with the experimental 
value of T = 14. 8^7 keV, within the uncertainty. Section III.E presents the details of the separation of the E-map 
into the main cluster and subcluster components. 

Section IV is dedicated to the convergence re-map from the weak and strong gravitational lensing survey of the Bullet 
Cluster 1E0657-558 from the November 15, 2006 data release re-map (Clowe et at, 2006b). Section IV. A presents a 
visualization of the re-map and some remarks on the evidence for dark matter or conversely, extra gravity. Section 
IV. B presents the MOG solution which uses a projective approximation for the density profile to facilitate numerical 
integration, although the full non-spherically symmetric expressions are provided in Section II. A. To a zeroth order 
approximation, we present the spherically symmetric solution, which does not fit the Bullet Cluster 1E0657-558. In 
our next approximation, the density profile of the Bullet Cluster 1E0657-558, p, is taken as the best-fit King /3-model 
(spherically symmetric), but the dimensionless gravitational coupling, Q, also assumed to be spherically symmetric, 
has a different center - in the direction closer to the subcluster. We determined that our best-fit re model corresponds to 
a location of the MOG center 140 kpc away from the main cluster E-map toward the sub-cluster E-map peak. Section 
IV. C presents the MOG prediction of the galaxy surface mass density, computed by taking the difference between the 
re-map data and our re-model of the E-map (ICM gas) data. Section IV. D presents the mass profile of dark matter 
computed by taking the difference between the re-map data and scaled E-map (ICM gas) data. This corresponds to 
the amount of dark matter (not a falsifiable prediction) necessary to explain the Bullet Cluster 1E0657-558 data using 
Einstein/Newton gravity theory. 

II. THE THEORY 

A. Modified Gravity (MOG) 

Analysis of the recent X-ray data from the Bullet Cluster 1E0657-558 (Markevitch, 2006), probed by our computa- 
tion in Modified Gravity (or MOG), provides direct evidence that the convergence re-map reconstructed from strong 
and weak gravitational lensing observations (Bradac et al., 2006; Clowe et al., 2006a, c) correctly accounts for all of the 
baryons in each of the main and subclusters including all of the galaxies in the regions near the main central dominant 
(cD) and the subcluster 's brightest central galaxy (BCG). The available baryonic mass, in addition to a second-rank 
skew symmetric tensor field (in MSTG), or massive vector field (in STVG), are the only properties of the system which 
contribute to the running gravitational coupling, G(r). It is precisely this effect which allows MOG to fulfill its require- 
ment as a relativistic theory of gravitation to correctly describe astrophysical phenomena without the necessity of dark 
matter (Moffat, 2005, 2006a,b). MOG contains a running gravitational coupling - in the infrared (IR) at astrophysical 
scales - which has successfully been applied to galaxy rotation curves (Brownstein and Moffat, 2006b), X-ray cluster 
masses (Brownstein and Moffat, 2006a), and is within limits set by solar system observations (Brownstein and Moffat, 
2006c). 

The weak field, point particle spherically symmetric acceleration law in MOG is obtained from the action principle 
for the relativistic equations of motion of a test particle in Moffat (2005, 2006b). The weak field point particle 
gravitational potential for a static spherically symmetric system consists of two parts: 

$(r) = $ N (r) + * y (r), (1) 

where 



(2) 
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and 



$y(r) = a 



exp(—fj,r) 



(3) 



denote the Newtonian and Yukawa potentials, respectively. M denotes the total constant mass of the system and /j, 
denotes the effective mass of the vector particle in STVG. The Poisson equations for $jv(r) and <&y(r) are given by 



V 2 $Ar(r) = -GooP(r), 



and 



F7 J - H 4 )*y(r) = — p(r), 



(4) 



(•5) 



respectively. For sufficiently weak fields, we can assume that the Poisson Equations (4) and (5) are uncoupled and 
determine the potentials <j>jv(r) and $y(r) for non-spherically symmetric systems, which can be solved analytically 
and numerically. The Green's function for the Yukawa Poisson equation is given by 



(V 2 - M 2 )Ay(r) = -<5 3 (r). 
The full solutions to the potentials are given by 



® N (v) = -G 00 J d 3 r' 



and 



$y(r) 



/ d 3 rV(r') 



, 3 ,exp(-/i|r- r'l)p(r') 



r — r' 



For a delta function source density 



p(r) = MS 3 (r), 



we obtain the point particle solutions of Equations (2) and (3). 
The modified acceleration law is obtained from 



a(r) = -V$ = -(V$ w (r) + V$ y (r)). 



Let us set 



Goo — Gn 



M 



Jd 3 r'p(r') 



1/2 



(6) 



(7) 



(8) 



(9) 



(10) 



(11) 



a = G N [M J dV P (r')] 1/2 , 



(12) 



where Mo is a constant and Gn denotes Newton's gravitational constant. From Equations (7), (8) and (10), we obtain 

a(r) = -G N [ d ^-^ 



r — r 



'1 3 



X<^ 1 + 



Mo 



1/2 



We can write this equation in the form: 



a(r) 



fd 3 r'p{r>) 
l-exp(- M |r-r'|)(l+ M |r-r'|) 



^^fG(r-r % 



(13) 



(14) 
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where 



G(r-r') = G N fl 



M 



1/2 



J d 3 r'p(r' 
x l-exp(-/i|r-r'|)(l + ^|r-r'|) 



(15) 



For a static spherically symmetric point particle system, we obtain using Equation (9) the effective modified 
acceleration law: 



a(r) 



G(r)M(r) 



(16) 



where 



am G,\l + J^ 



1 - exp(-r/r ) H 



(17) 



Here, M is the total baryonic mass of the system and we have set fi = 1/ro and r is a distance range parameter. We 
observe that G(r) — > G^r as r — > 0. 
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FIG. 2: 

Running Scales in Modified Gravity. 

A plot of the MOG mass scale, 
Mo, vs. the MOG range pa- 
rameter, ro compiled from 
published fits to galaxy rota- 
tion curves and X-ray cluster 
masses (Brownstein and Moffat, 
2006a,b) The values of M vs.ro 
are plotted for galaxies with green 
squares, dwarf galaxies with black 
circles, clusters of galaxies with 
blue crosses, and dwarf clusters 
of galaxies with red triangles. 



For the spherically symmetric static solution in MOG, the modified acceleration law Equation (16) is determined by 
the baryon density and the parameters Gn, Mq and ro. However, the parameters Mo and ro are scaling parameters 
that vary with distance, r, according to the field equations for the scalar fields u>(r) oc Mo(r) and fi(r) = l/r (r), 
obtained from the action principle (Moffat, 2006b). In principle, solutions of the effective field equations for the 
variations of u>{r) and /j,(r) can be derived given the potentials V(u>) and V(/j,) in Equations (24) and (26) of Moffat 
(2006b) However, at present the variations of w(r) and fi(r) are empirically determined from the galaxy rotation curve 
and X-ray cluster mass data (Brownstein and Moffat, 2006a,b). 

In Figure 2, we display the values of Mo and ro for spherically symmetric systems obtained from the published 
fits to the galaxy rotation curves for dwarf galaxies, elliptical galaxies and spiral galaxies, and the fits to (normal 
and dwarf) X -ray cluster masses. A complete continuous relation between Mo and ro at all mass scales needs to be 
determined from the MOG that fits the empirical mass scales and distance scales shown in Figure 2. This will be a 
subject of future investigations. 

The modifications to gravity leading to Equation (17) would be negated by the vanishing of either Mo — > or the 
r — ► oo infrared limit. These scale parameters are not to be taken as universal constants, but are source dependent 
and scale according to the system. This is precisely the opposite case in MOND, where the Milgrom acceleration, 
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aoMiigiom = 1-2 x 10 -8 cm s~ 2 (Milgrom, 1983), is a phenomenologically derived universal constant - arising from 
a classical modification to the Newtonian potential. Additionally, MOND has an arbitrary function, fi(x), which 
should be counted among the degrees of freedom of that theory. Conversely, MOG does not arise from a classical 
modification, but from the equations of motion of a relativistic modification to general relativity. 

The cases we have examined until now have been modeled assuming spherical symmetry. These include 101 
galaxy rotation curves (Brownstein and Moffat, 2006b), 106 X-ray cluster masses (Brownstein and Moffat, 2006a), 
and the gravitational solution to the Pioneer 10/11 anomaly in the solar system (Brownstein and Moffat, 2006c). In 
applications of MOG, we vary the gravitational coupling G, the vector field </> M coupling to matter and the effective mass 
\x of the vector field according to a renormalization group (RG) flow description of quantum gravity theory formulated 
in terms of an effective classical action (Moffat, 2005, 2006b). Large infrared renormalization effects may cause the 
effective coupling constants to run with momentum. A cutoff leads to spatially varying values of G, M and r and 
these values increase in size according to the mass scale and distance scale of a physical system (Rcutcr and Weyer, 
2006). 

The spatially varying dimcnsionless gravitational coupling is given by 



where 



• G N = 6.6742 x 10- 11 m 3 /kg s 2 is the ordinary (terrestrial) Newtonian gravitational constant measured 
experimentally 2 , 



M(r) is the total (ordinary) mass enclosed in a sphere or radius, r. This may include all of the visible (X-ray) 
ICM gas and all of the galactic (baryonic) matter, but none of the non-baryonic dark matter. 



Mo is the MOG mass scale (usually measured in units of [M©]), 
• ro is the MOG range parameter (usually measured in units of [kpc]). 
The dimcnsionless gravitational coupling, Q(r), of Equation (18), approaches an asymptotic value as r — > oo: 

g °° =i+ \/W' (i9) 

where M is the total baryonic mass of the system. Conversely, taken in the limit of r <C r , G(r) — > Gn or Q(r) — ► 1 
down to the Planck length. 

In order to apply the MOG model of Equation (18) to the Bullet Cluster 1E0657-558, we must first generalize the 
spherically symmetric case. Our approach follows a sequence of approximations: 

i. Treat the subcluster as a perturbation, and neglect it as a zeroth order approximation. 

ii. Treat the subcluster as a perturbation, and shift the origin of the gravitational coupling, Q{r) toward the 
subcluster (toward the center-of-mass of the system). 

iii. Use the concentric cylinder mass M{R) as an approximation for M(r), and shift that toward the subcluster 
(toward the MOG center). 

iv. Treat the subcluster as a perturbation, and utilize the isothermal sphere model to approximate M(R) and shift 
that toward the subcluster (toward the center-of-mass of the system - where the origin of Q(r) is located). 



B. The King /3-Model for the S-map 

Starting with the zeroth order approximation, we neglect the subcluster, and perform a best-fit to determine a 
spherically symmetric King /3-model density of the main cluster. We assume the main cluster gas (neglecting the 
subcluster) is in nearly hydrostatic equilibrium with the gravitational potential of the galaxy cluster. Within a few 



2 NIST 2002 CODATA value. 
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core radii, the distribution of gas within a galaxy cluster may be fit by a King "/3-model" . The observed surface 
brightness of the X-ray cluster can be fit to a radial distribution profile (Chandrasckhar, 1960; King, 1966): 



I(r) = Jo 



1+1 



-i -3/3+1/2 



(20) 



resulting in best-fit parameters, (3 and r c . A deprojection of the /3-model of Equation (20) assuming a nearly isothermal 
gas sphere then results in a physical gas density distribution (Cavaliere and Fusco-Femiano, 1976): 



P(r) = Pq 



, r 
1+ - 



-3/3/2 



(21) 



where p(r) is the ICM mass density profile, and po denotes the central density. The mass profile associated with this 
density is given by 

pr 

(22) 



M(r) = 4tt / p(r')r' dr', 
Jo 

where M(r) is the total mass contained within a sphere of radius r. Galaxy clusters are observed to have finite spatial 
extent. This allows an approximate determination of the total mass of the galaxy cluster by first solving equation 
(21) for the position, r out , at which the density, p(r ou t), drops to « 10~ 28 g/cm 3 , or 250 times the mean cosmological 
density of baryons: 



r out — r c 



Po 



2/3/3 



10~ 28 g/crn 



1/2 



(23) 



Then, the total mass of the ICM gas may be taken as M gas ss M(r out ): 



M gas = 4tt 











1 + 




/ Po 




h 







(24) 



To make contact with the experimental data, we must calculate the surface mass density by integrating p(r) of 
Equation (21) along the line-of-sight: 



p(x,y,z)dz, 



where 



^out 



x — y 



2 , 2 , 2 1 -3/3/2 

x + y + z A 



Substituting Equation (21) into Equation (25), we obtain 



= po 



This integral becomes tractable by making a substitution of variables: 

u =1-1 ^ — , 



dz. 



(25) 



(26) 



(27) 



(28) 



so that 



^(x,y) = Po 



Po 

U 3P 



2out 

"2out 
z out 



z 
ur. 



-3/3/2 



dz 

-3/3/2 



,,3/3 Zmtr 



1 3 

2' 2 







dz 

( ^out 

V ur c 



(29) 



10 



where we have made use of the hypergeometric function, F([a,b], [c],z). Substituting Equation (28) into Equation 
(29) gives 



T,(x,y) = 2p z out 1 + 



x 2 + y 2 " 



N "3/3/2 / 






"3" 


) F { 






2 



x 2 + y 2 + r 2 



We next define 

So = S(0, 0) = 2p z out F 
which we substitute into Equation (30), yielding 

E(x,y) = S ( 1 




(30) 



(31) 



x 2 +y 2 Y W/2F {[hMi 



21 ' x 2 +y^+r^ 



(32) 



In the limit z out S> r c , the Hypergeometric functions simplify to Y functions, and Equations (31) and (32) result in 
the simple, approximate solutions: 



S = Vnpor c 



-(^) 

r(§/3) 



and 



(33) 



(34) 



T,(x,y) = S I 1 

which we may, in principle, fit to the S-map data to determine the King /3-model parameters, {3, r c and po- 

C. Deriving the Weighted Surface Mass Density from the Convergence s-map 

The goal of the strong and weak lensing survey of Bradac et al. (2006); Clowe et al. (2006a,c) was to obtain a 
convergence K-map by measuring the distortion of images of background galaxies (sources) caused by the deflection 
of light as it passes the Bullet Cluster 1E0657-558 (lens). The distortions in image ellipticity are only measurable 
statistically with large numbers of sources. The data were first corrected for smearing by the point spread function 
in the image, resulting in a noisy, but direct, measurement of the reduced shear g = 7/(1 — n). The shear, 7, is 
the anisotropic stretching of the galaxy image, and the convergence, k, is the shape-independent change in the size 
of the image. By recovering the K-map from the measured reduced shear field, a measure of the local curvature is 
obtained. In Einstein's general relativity, the local curvature is related to the distribution of mass/energy, as it is in 
MOG. In Newtonian gravity theory the relationship between the K-map and the surface mass density becomes very 
simple, allowing one to refer to n as the scaled surface mass density (see, for example, Chapter 4 of Peacock (2003) 
for a derivation): 



n(x,y) 



where 



is the surface mass density, and 



&■ D s S c 



^(x,y) = / p(x,y,z)dz, 



4nG N A As 



3.1 x 10 9 M /kpc 



(35) 



(36) 



(37) 
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is the Newtonian critical surface mass density (with vanishing shear), D s is the angular distance to a source (back- 
ground) galaxy, D\ is the angular distance to the lens (Bullet Cluster 1E0657-558), and D\ s is the angular distance 
from the Bullet Cluster 1E0657-558 to a source galaxy. The result of Equation (37) is equivalent to 

« 540 kpc. (38) 

Since there is a multitude of source galaxies, these distances become "effective" , as is the numeric value presented in 
Equation (37), quoted from (Clowe et al., 2004) without estimate of the uncertainty. In fact, due to the multitude of 
sources in the lensing survey, both D s and D\ s are distributions in (x,y). However, it is common practice to move 
them outside the integral, as a necessary approximation. 

We may obtain a similar result in MOG, as was shown in Moffat (2006a), by promoting the Newtonian gravitational 
constant, Gjv to the running gravitational coupling, G(r), but approximating G(r) as sufficiently slow- varying to allow 
it to be removed from the integral. We have in general, 

«*>V) = J ^^p(*,y,*)d Z = f^, (39) 

and we assume that 

c 2 D s Gm E c 
47rG(r) D\D\ S G(r) y(r) 

where we applied Equation (18). However, Equations (39) and (40) are only valid in the thin lens approximation. For 
the Bullet Cluster 1E0657-558, the thin lens approximation is inappropriate, and we must use the correct relationship 
between k and E: 

= / ^^K*,v,*)*< - ^ (4i) 

where 

t(x,y) = J g(r)p(x,y,z)dz, (42) 

is the weighted surface mass density (weighted by the dimensionless gravitational coupling Q(r) of Equation (18)), 
and E c is the usual Newtonian critical surface mass density Equation (37). 

For the remainder of this paper, we will use Equations (18), (37), (41) and (42) to reconcile the experimental 
observations of the gravitational lensing K-map with the X-ray imaging E-map. We can already see from these 
equations, how in MOG the convergence, n, is now related to the weighted surface mass density, E, so 

... k is no longer a measurement of the surface density, but is a non-local function whose overall level 
is still tied to the amount of mass. For complicated system geometries, such as a merging cluster, the 
multiple peaks can deflect, suppress, or enhance some of the peaks (Clowe et at, 2006c). 



III. THE SURFACE DENSITY MAP FROM X-RAY IMAGE OBSERVATIONS 

A. The E-map 

With an advance of the November 15, 2006 data release (Clowe et al., 2006b), we began a precision analysis to 
model the gross features of the surface density E-map data in order to gain insight into the three-dimensional matter 
distribution, p(r), and to separate the components into a model representing the main cluster and the subcluster - 
the remainder after subtraction. 

The E-map is shown in false colour in Figure 3. There are two distinct peaks in the surface density E-map - the 
primary peak centered at the main cluster, and the secondary peak centered at the subcluster. The main cluster gas 
is the brightly glowing (yellow) region to the left of the subcluster gas, which is the nearly equally bright Shockwave 
region (arrowhead shape to the right). The K-map observed peaks, the central dominant (cD) galaxy of the main 
cluster, the brightest cluster galaxy (BCG) of the subcluster, and the MOG predicted gravitational center are shown 
in Figure 3 for comparison. J2000 and map (x,y) coordinates are listed in Table 1. 

We may now proceed to calculate the best-fit parameters, (3, r c and po of the King /3-model of Equations (33) 
and (34), by applying a nonlinear least-squares fitting routine (including estimated errors) to the entire E-map, or 
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FIG. 3: Surface density E-map. 

Data reconstructed from X-ray imaging observations of the Bullet Cluster 1E0657-558, November 15, 2006 data 
release (Clowe et al., 2006b). E-map observed peaks (local maxima) and K-map observed peaks are shown for 
comparison. The central dominant (cD) galaxy of the main cluster, the brightest cluster galaxy (BCG) of 
the subcluster, and the MOG predicted gravitational center are shown. The colourscale is shown at the bot- 
tom, in units of 10 15 MQ/pixel 2 . The resolution of the E-map is 8.5 kpc/pixel, based upon the mea- 
sured redshift distance ~ 260.0 kpc/arcminute of the Bullet Cluster 1E0657-558 (Bradac et al., 2006). The 
scale in kpc is shown at the top of the figure. J2000 and map (x,y) coordinates are listed in Table 1. 



alternatively by performing the fit to a subset of the E-map on a straight-line connecting the main cluster S-map 
peak (R = 0) to the main cD, and then extrapolating the fit to the entire map. This reduces the complexity of the 
calculation to a simple algorithm, but is not guaranteed to yield a global best-fit. However, our approximate best-fit 
will prove to agree with the S-map everywhere, except at the subcluster (which is neglected for the best-fit). 

The scaled surface density E-map data is shown in solid red in Figure 4. The unmodeled peak (at R ~ 300 kpc) is 
due to the subcluster. The best-fit to the King /3-model of Equation (34) is shown in Figure 4 in short-dashed blue, 
and corresponds to 

13 = 0.803 ±0.013, (43) 
r c = 278.0 ±6.8 kpc, (44) 

where the value of the E-map at the main cluster peak is constrained to the observed value, 

E = 1.6859 x 10 10 ^% f V = 2-3181 x lO 8 M /kpc 2 . (45) 

pixel V 8 - 528 k P c / 

We may now solve Equation (33) for the central density of the main cluster, 

v y ( -8) 

po = , \ = 3.34 x 10 5 M /kpc 2 . (46) 

V7rr c r " ' 1 



(2^1) 
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Observation 




J2000 Coordinates 




E-map 


ft-map 






RA 




Dec 




{x, y) 


(x, y) 


Main cluster E-map peak 


06 


58 


31.1 


-55 


56 


53.6 


(89, 89) 


(340, 321) 


Subcluster E-map peak 


06 


58 


20.4 


-55 


56 


35.9 


(135, 98) 


(365, 326) 


Main cluster K-map peak 


06 


58 


35.6 


-55 


57 


10.8 


(70, 80) 


(329, 317) 


Subcluster K-map peak 


06 


58 


17.0 


-55 


56 


27.6 


(149, 102) 


(374, 327) 


Main cluster cD 


06 


58 


35.3 


-55 


56 


56.3 


(71,88) 


(330, 320) 


Subcluster BCG 


06 


58 


16.0 


-55 


56 


35.1 


(154, 98) 


(375, 326) 


MOG Center 


06 


58 


27.6 


-55 


56 


49.4 


(105,92) 


(348, 322) 



TABLE 1: The J2000 sky coordinates of the Bullet Cluster 1E0657-558, November 15, 2006 data release (Clowe er a/., 2006b): 
Main cluster and subcluster E-map and K-map observed peaks, the central dominant (cD) galaxy of the 
main cluster, the brightest cluster galaxy (BCG) of the subcluster, and the MOG predicted gravitational cen- 
ter. The resolution of the E-map and K-map are 8.5 kpc/pixel, and 15.4 kpc/pixel, respectively, based upon 
the measured redshift distance ~ 260.0 kpc/arcminute of the Bullet Cluster 1E0657-558 (Bradac et ai, 2006). 



The values of the parameters, (3, r c and po completely determine the King /3-model for the density, p(r), of Equation 
(21) of the main cluster X-ray gas. We provide a comparison of the full E-map data (Figure 5a, in red) with the 
result of the surface density S-map derived from the best-fit King /3-model to the main cluster (Figure 5b, in blue) . 
Substituting Equations (43), (44) and (46) into Equation (23), we obtain the main cluster outer radial extent, 

r out = 2620 kpc, (47) 

the distance at which the density, p{r out ), drops to sa 10~ 28 g/cm 3 , or 250 times the mean cosmological density of 
baryons. The total mass of the main cluster may be calculated by substituting Equations (43), (44) and (46) into 
Equation (24): 

M gas = 3.87 x 10 14 M Q . (48) 



0.08 



0.06 - 




0.01 



0.02 



1000 



FIG. 4: King /3-model fit to scaled E-map. 

A cross-section of the scaled surface density E- 
map data reconstructed from X-ray imaging ob- 
servations of the Bullet Cluster 1E0657-558, on 
a straight-line connecting the main X-ray cluster 
peak (R = kpc) to the main central domi- 
nant (cD) galaxy (R w -150 kpc). The Main 
cluster E-map data, taken from the November 15, 
2006 data release (Clowe et ai, 2006b), is shown 
in solid red, and the surface density E-map ac- 
cording to the best-fit King /3-model (neglecting 
the subcluster) of Equation (34) is shown in short- 
dashed blue. The unmodeled peak (at 7? ~ 
300 kpc) is due to the subcluster. We used the 
Clowe et al. (2004) value for the Newtonian crit- 
ical surface mass density (with vanishing shear), 
E c = 3.1 x 10 9 M Q /kpc 2 . J2000 and map (x,y) 
coordinates are listed in Table 1. The best-fitting 
King /3-model parameters are listed in Table 2. 
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120 ieo 200 



(a) E-map 



(b) King /3-model 



FIG. 5: Comparison of the E-map data and the best-fit King /3-model. 

The surface density E-map data reconstructed from X-ray imaging observations of the Bullet Clus- 
ter 1E0657-558, November 15, 2006 data release (Clowe el ai, 2006b), is shown in red (a). The 
surface density E-map according to the best-fit King /3-model (neglecting the subcluster) of Equa- 
tion (34) is shown in blue (b). The best-fitting King /3-model parameters are listed in Table 2. 



B. The Gravitational Coupling for the Main cluster 

As discussed in Section II. A, in order to apply the MOG model of Equation (18) to the Bullet Cluster 1E0657-558, 
we must first generalize the spherically symmetric case by treating the subcluster as a perturbation. In the zeroth 
order approximation, we begin by neglecting the subcluster. According to the large (> 100) galaxy cluster survey 
of Brownstein and Moffat (2006a), the MOG mass scale, Mq, is determined by a power-law, depending only on the 
computed total cluster mass, M gas : 

/ \ 0.39±0.10 

Mo = (60.4 ± 4.1) x 10 14 M [jq^) ■ (49) 

Substituting the result of Equation (48) for the main cluster into Equation (49), we obtain 

M a = 1.02 x 10 16 M @ , (50) 

and substituting the result of Equation (48) for M and Equation (50) for the main cluster of Bullet Cluster 1E0657-558 
into Equation (19), we obtain 

£oo = 6.14. (51) 



King Model 


Main Cluster 


Subcluster 


P 


0.803 ±0.013 




r c 


278.0 ±6.8 kpc 




PQ 


3.34 x 10 5 M e /kpc 2 




Afgas 


3.87 x 10 14 M Q 


2.58 x 10 13 Mq 


^out 


2620 kpc 





TABLE 2: The best fit King /3-model parameters. 

Using a nonlinear least-squares fitting routine, 
Equations (33) and (34) are fit to the E-map 
data for the main cluster (neglecting the sub- 
cluster). Results of our calculation for the 
main and subcluster mass, Mg as , and cluster 
outer radial extent, r out are listed below. The 
best-fit is shown in Figures 4 and 5b in blue, 
and may be compared to the data, in red. 
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MOG 


Main Cluster 


Subcluster 


Best-fit to K-map 


M 
ro 


1.02 x 10 16 M & 
139.2 kpc 


3.56 x 10 15 M Q 


3.07 x 10 15 M Q 
208.8 kpc 


Goo 

T 


6.14 
15.5 ±3.9 keV 




3.82 



FIG. 6: The MOG mass scale, Mo, vs. M gas . 
A graphical calculation based on a power-law 
running of the MOG mass scale, Mo, depend- 
ing on the total ICM gas mass, M gas , de- 
termined by the large (> 100) galaxy cluster 
survey of Brownstein and Moffat (2006a). The 
computed values of M gas and Mo for the 
main and subcluster are listed in Table 3. 



TABLE 3: MOG parameters. 

Results of calculation of the power-law running of the MOG 
mass scale, Mo, depending on the computed total cluster 
mass, Mg as , according to the large (> 100) galaxy cluster 
survey of Brownstein and Moffat (2006a) for the main and 
subcluster (shown graphically in Figure 6), and the overall 
best-fit to the K-map; and the results for the MOG range 
parameter, ro. Results of the calculations for the asymp- 
totic value of the dimensionless gravitational coupling, 
Goo, and the calculated X-ray temperature of the unper- 
turbed isothermal sphere for the main cluster, T, are listed. 



From Brownstein and Moffat (2006a), the MOG range parameter, ro, depends only on the computed outer radial 
extent, r out : 



from which we obtain 



= rout/10, r out < 650 kpc, 
r = 139.2 kpc, r out > 650 kpc, 

r = 139.2 kpc, 



(52) 



(53) 



for the main cluster. 

The best-fitting King /3-model parameters for the main cluster are listed in Table 2, and the MOG parameters 
for the main cluster are listed in Table 3. A plot of the dimensionless gravitational coupling of Equation (18) using 
the MOG parameter results of Equations (50) and (53) for the main cluster (neglecting the subcluster) is plotted in 
Figure 7a using a linear scale for the r-axis, and in Figure 7b on a logarithmic scale for the r-axis. 



C. The Cylindrical Mass Profile 

y) is an integrated density along the line-of-sight, z, as shown in Equation (36). By summing the S-map pixel- 
by-pixel, starting from the center of the main cluster E-map peak, one is performing an integration of the surface 
density, yielding the mass, 



M(R) 



Y<{R')R'dR', 



(54) 
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1000 1500 

T [kpc] 

(a) Linear Scale, r-axis. 




(b) Logarithmic Scale, r-axis. 



FIG. 7: Plot of the dimensionless gravitational coupling, Q(r) = G(t)/Gn, of Equation (18) vs. the distance, r, in kpc. 

Shown in linear scale for the r-axis, (a), and in logarithmic scale, (b). J2000 and map (x,y) coordinates of the MOG 
center (R = kpc) are listed in Table 1, and located with respect to the E-map and K-map in Figures 3 and 12, respec- 
tively. The running gravitational coupling, G(r), shown here, corresponds to our best-fit to the K-map, listed in Table 3. 



sr io 12 




100 
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1000 



FIG. 8: Cylindrical mass profile 

The surface density E-map reconstructed from X- 
ray image observations of the Bullet Cluster 1E0657- 
558, November 15, 2006 data release (Clowe et ai, 
2006b) is integrated about the main cluster E-map 
peak (R = 0), shown in solid red. The cylindri- 
cal mass profile derived from the best-fit King /3- 
model is shown in short-dashed blue, demonstrating 
a good agreement with the data. The long-dash ma- 
genta line is the result of the galaxy subtraction of 
Equations (78) and (79). The lack of galactic mass 
in the vicinity of the main cluster E-map peak at 
R — is due to the merger. The dash-dot black 
line is the result of the dark matter subtraction of 
Equation (83). The dark matter dominates the ICM 
cluster gas by a factor > 3. J2000 and map (x,y) 
coordinates are listed in Table 1. The best-fitting 
King /3-model parameters are listed in Table 2. 



enclosed by concentric cylinders of radius, R = \J x 1 + y 2 . We performed such a sum over the E-map data, and 
compared the result to the best-fit King /3-model derived in Section III. A, with the parameters listed in Table 2. The 
results are shown in Figure 8. The fact that the data and model are in good agreement provides evidence that the King 
/3-model is valid in all directions from the main cluster S-map peak, and not just on the straight line connecting the 
peak to the main cD, where the fit was performed. The model deviates slightly from the data, underpredicting M(R) 
for R > 350 kpc, which may be explained by the subcluster (which is included in the data, but not the model). The 
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integrated mass profile arising from the dark matter analysis of the K-map of Equation (83) is shown for comparison 
on the same figure. 

The ratio of dark matter to ICM gas for the main cluster is > 3, which is significantly less than the average 
cosmological ratio of 5.68 discussed in Section LA, whereas one should expect the highest mass to light ratios from 
large, hot, luminous galaxy clusters, and the Bullet Cluster 1E0657-558 is certainly one of the largest and hottest. 



D. The Isothermal Spherical Mass Profile 



For a spherical system in hydrostatic equilibrium, the structure equation can be derived from the collisionless 
Boltzmann equation 



O(rK) + — — (cr r - o- M ) = -p{r)- 



dr 



dr 



(55) 



where <£>(r) is the gravitational potential for a point source, ay and (Je,<j> are mass- weighted velocity dispersions in the 
radial (r) and tangential {9, <f)) directions, respectively. For an isotropic system, 



c r — ve,c/>- 

The pressure profile, P(r), can be related to these quantities by 

P(r) = <j 2 r p(r). 

Combining equations (55), (56) and (57), the result for the isotropic sphere is 



dP(r) 



-p(r) 



d$(r) 



dr dr 

For a gas sphere with temperature profile, T(r), the velocity dispersion becomes 

2 kT(r) 
a r = , 

PATflp 



(56) 
(57) 

(58) 
(59) 



where k is Boltzmann's constant, pa ~ 0.609 is the mean atomic weight and m p is the proton mass. We may now 
substitute equations (57) and (59) into equation (58) to obtain 



d (kT{r) 



dr \pAm v 



P( r ) = ~p{r) 



d<S>(r) 
dr 



(60) 



Performing the differentiation on the left-hand side of equation (58), we may solve for the gravitational acceleration: 



_ d#(r) 

a(r) = dV 

kT(r) 
p A m p r 



d\n(p(r)) d\n(T(r)) 



dln(r) 



dln(r) 



(61) 



For the isothermal isotropic gas sphere, the temperature derivative on the right-hand side of equation (61) vanishes 
and the remaining derivative can be evaluated using the /3-model of equation (21): 



a(r) = - 



3/3/cT 



PArrip \r z + r 



(62) 



The dynamical mass in Newton's theory of gravitation can be obtained as a function of radial position by replacing 
the gravitational acceleration with Newton's Law: 



so that equation (61) can be rewritten as 

M N (r) = 



a N (r) 



r kT 
G N p A m v 



G N M(r) 



dhi(p(r)) dln{T(r)) 



dhi(r) dln(r) 



(63) 



(64) 
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and the isothermal /3-model result of equation (62) can be rewritten as 



M N (r) 



3(3kT 
[i A m v G N \r 2 + r 



(65) 



Similarly, the dynamical mass in MOG can be obtained as a function of radial position by substituting the MOG 
gravitational acceleration law (Brownstcin and Moffat, 2006a,b; Moffat, 2005, 2006b) 



a(r) = - 



G(r)M 



so that our result for the isothermal /3-model becomes 



M MO G{r) 



3(3kT 



HArripGir) \r 2 + r 2 

We can express this result as a scaled version of equation (64) or the isothermal case of equation (65): 

M N (r) 



(66) 



(67) 



A4iog(»") 




M 



Mmog(?") 



1 - exp(-r/r ) 1 + 



ro 



M N (r), 



(68) 



where we have substituted Equation (18) for Q{r). Equation (68) may be solved explicitly for MmogW by squaring 
both sides and determing the positive root of the quadratic equation. 

The scaling of the Newtonian dynamical mass by Q (r) according to Equation (68) solved the dark matter problem 
for the galaxy clusters of the > 100 galaxy cluster survey of Brownstein and Moffat (2006a). The unperturbed Bullet 
Cluster 1E0657-558 is no exception! In Figure 9, we plotted the MOG and the Newtonian dynamical masses, -Mmog(^ ) 
and M^fr), respectively, and compared it to the spherically integrated best- fit King /3-model for the main cluster gas 
mass of Equations (21) and (22) using the parameters listed in Table 2. The MOG temperature prediction due to the 
best-fit is listed in Table 4 and compared to experimental values. 

As shown in Figure 9, across the range of the r-axis, and throughout the radial extent of the Bullet Cluster 1E0657- 
558, the la correlation between the gas mass, M(r) and the MOG dynamical mass, MmogMj provides excellent 













































,-'* ■' 

■ * ■ / / , 








• ,' 
















//'/ 








/ ■' ' / 








' '' *jr- ■'' 
















'/ 








>/ ' 








■ fs ■ ' 








■ ks/ 








' V ' 






















*'* -' 








, * ■' // ■ 








, ' ■ 








/ 1- 








'/ f 








— '/ ' 








. ■/ / 








■/ / 
























-V 
V 








^ 1 1 









10 



1(10 
[kpc] 



1000 



FIG. 9: The spherically 
integrated mass profile of the main cluster ICM gas. 

The mass profile integrated King /3-model, M gas (r), 
in from the S-map peak, is shown in long-dashed 
red. derived from the surface density E-map data re- 
constructed from X-ray imaging observations of the 
Bullet Cluster 1E0657-558, November 15, 2006 data 
release (Clowe et at, 2006b). The surface density 
E-map according to the best-fit King /3-model (ne- 
glecting the subcluster) of Equation (34) is shown in 
black (with uncertainty). The Newtonian dynami- 
cal mass (the dark matter paradigm) is shown in 
short-dashed blue. All masses are presented in units 
of Mq, and distances, r, in kpc. The strong cor- 
respondence between the data (red) and the MOG 
fit (black) across distance scales in the cluster from 
tens to thousands of kpc. is significant. This 
is the first evidence that the X-ray imaging data 
for the cluster is behaving similar to the large (> 
100) galaxy cluster survey of Brownstein and Moffat 
(2006a). J2000 and map (x,y) coordinates are 
listed in Table 1. The best-fitting King /3-model 
parameters are listed in Table 2. The best- 
fitting MOG parameters are listed in Table 3. 
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Year 


Source - Theory or Experiment 


1 ^Kev ) 


% error 




MOG Prediction 


15.5 ±3.9 




2002 


accepted experimental value 


UAt\l 


4.5 


1999 


ASCA+ROSAT fit 


14.5+ 2 7 ° 


6.5 


1998 


ASCA fit 


17.4 ± 2.5 


12.3 



TABLE 4: The MOG best-fit tempera- 
ture is consistent with the experimental values for the main 
cluster isothermal temperature, the 1999 ASCA+ROSAT 
fit, and the 1998 ASCA fit (Markevitch et al, 2002). 



agreement between theory and experiment. The same cannot be said of any theory of dark matter in which the 
X-ray surface density map is negligible in relation to the DM density. So dark matter makes no prediction for the 
isothermal temperature, which has been measured to reasonable precision for many clusters, but simply "accounts for 
missing mass." Since there is no mysterious missing mass in MOG, the prediction should be taken seriously as direct 
confirmation of the theory, or at least as hard evidence for MOG. 



E. The Subcluster Subtraction 
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Main cluster E-map peak 
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Main cluster K-map peak 



MOG Center 




2E-06 4E-06 6E-06 8E-06 1 E-05 1 .2E-05 1 .4E-05 1 .6E-05 



FIG. 10: The subcluster subtracted surface density E-map 

Main cluster and subcluster n-map observed peaks (local maxima) and E-map observed peaks are shown 
for comparison. The central dominant (cD) galaxy of the main cluster, the brightest cluster galaxy 
(BCG) of the subcluster, and the MOG predicted gravitational center are shown. The colourscale is 
set to agree with Figure 3 for comparison. J2000 and map (x,y) coordinates are listed in Table 1. 



Provided the surface density E-map derived from the best-fit King /3-model of Equation (34) is sufficiently close 
to the E-map data (consider Figures 4 and 5) then the difference between the data and the /3-model is the surface 
density E-map due to the subcluster. Our subcluster subtraction, shown in Figures 10 and 11a, is based upon a 
high precision (x 2 < 0.2) best-fit King /3-model to the main cluster, which agrees with the main cluster surface mass 
E-map (data) within 1% everywhere and to a mean uncertainy of 0.8%. The subcluster subtraction is accurate down 
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to p — 10 -28 g/cm 3 ~ 563.2 M /pc 3 baryonic background density. After subtraction, the subcluster S-map peak 
takes a value of 1.30 x 10 8 M /kpc 2 , whereas the full S-map has a value of 2.32 x 10 8 M /kpc 2 at the subcluster 
S-map peak. Thus the subcluster (at its most dense position) provides only w 56% of the X-ray ICM, the rest is 
due to the extended distribution of the main cluster. The subcluster subtraction surface density S-map shown in 
Figure 10 uses the same colorscale as the full S-map shown in Figure 3, for comparison. Figure 11 is a stereogram 
of the subcluster subtracted surface density S-map and the subcluster superposed onto the surface density S-map of 
the best-fit King /3-model to the main cluster. There is an odd X-ray bulge in Figures 10 and lib which may be an 
artifact of the subtraction, or perhaps evidence of an as yet unidentified component. 

Since the outer radial extent of the subluster gas is less than 400 kpc, the S-map completely contains all of the 
subcluster gas mass. By summing the subcluster subtracted S-map pixel-by-pixel over the entire S-map peak, one is 
performing an integration of the surface density, yielding the total subcluster mass. We performed such a sum over 
the subcluster subtracted S-map data, obtaining 

M gas = 2.58 x 10 13 M , (69) 

for the mass of the subcluster gas, which is less than 6.7% of the mass of main cluster gas (the best-fitting King 
/3-model parameters are listed in Table 2.) This justifies our initial assumption that the subcluster may be treated as 
a perturbation in order to fit the main cluster to the King /3-model. Our subsequent analysis of the thermal profile 
confirms that the main cluster X-ray temperature is nearly isothermal, lending further support to the validity of the 
King /3-model. 




V 



(a) The subcluster subtracted surface density E-map (b) The subcluster subtracted E-map super- 

posed onto the best-fit King /3-model of the main clus- 
ter. Note the odd X-ray Bulge, here, and in Figure 10. 



FIG. 11: The subcluster subtracted surface density E-map. 

The blue surface represents the E-map due to the integrated (line-of-sight) King /3-model fit to main cluster (as shown in 
blue in Figure 5b). The green surface is the contribution to the E-map from the subcluster (as calculated from our model). 
J2000 and map (x,y) coordinates are listed in Table 1. The best-fitting King /3-model parameters are listed in Table 2. 

We may now calculate the MOG mass scale, M , of the subcluster by substituting the subcluster gas mass, M gas , 
of Equation (69) into the power-law relation of Equation (49), yielding, 

M = 3.56 x 10 15 M , (70) 

as shown in Figure 6. The computed values of M gas and Mq for the main and subcluster are listed in Table 3. 



21 



IV. THE CONVERGENCE MAP FROM LENSING ANALYSIS 



A. The re-Map 




FIG. 12: The surface density ft-map reconstructed from strong and weak gravitational lensing. 

Main cluster and subcluster of the Bullet Cluster 1E0657-558, November 15, 2006 data release ft-map (Clowe et al., 
2006b) observed peaks (local maxima) and E-map observed peaks are shown for comparison. The cen- 
tral dominant (cD) galaxy of the main cluster, the brightest cluster galaxy (BCG) of the subcluster, and 
the MOG predicted gravitational center are shown. J2000 and map (x,y) coordinates are listed in Table 1. 



As tempting as it is to see the convergence K-map of Figure 12 - a false color image of the strong and weak 
gravitational lensing reconstruction (Bradac et al., 2006; Clowe et al., 2006a,c) - as a photograph of the "curvature" 
around the Bullet Cluster 1E0657-558, it is actually a reconstruction of all of the bending of light over the entire 
distance from the lensing event source toward the Hubble Space Telescope. The source of the K-map is oc J Gnp{t), 
along the line-of-site, as in Equation (35), oc G(r) J p(r) as in Equation (39), or tx J G(r)p(r) as in Equation (41). 
For the Bullet Cluster 1E0657-558, we are looking along a line-of-sight which is at least as long as indicated by a 
redshift of z — 0.296 (Gpc scale). The sources of the lensing events are in a large neighbourhood of redshifts, an 
estimated z — 0.85 ± 0.15. This fantastic scale (several Gpc) is naturally far in excess of the distance scales involved 
in the X-ray imaging surface density E-map. It is an accumulated effect, but only over the range of the X-ray source 
- as much as 2.2 Mpc. A comparison of these two scales indicates that the distance scales within the E-map are 10~ 3 
below the Gpc's scale of the K-map. 

Preliminary comments on the November 15, 2006 data release (Clowe et al., 2006b): 

• The conclusion of Bradac et al. (2006); Clowe et al. (2006a), that the K-map shows direct evidence for the 
existence of dark matter may be premature. Until dark matter has been detected in the lab, it remains an open 
question whether a modified gravity theory, such as MOG, can account for the K-map without nonbaryonic dark 
matter. MOG, due to the varying gravitational coupling, Equation (18), gives the Newtonian 1/r 2 gravitational 
force law a considerable boost - "extra gravity" as much as Goo ~ 6 for the Bullet Cluster 1E0657-558. 
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• It may be feasible that the mysterious plateau in the north-east corner of the K-map is from some distribution 
of mass unrelated to the Bullet Cluster 1E0657-558. This "background curvature" contribution to the K-map is 
one part of the uncertainty in the reconstruction. The second, dominant source of uncertainty certainly must be 
the estimate of the angular diameter distances between the source of the lcnsing event and the Bullet Cluster 
1E0657-558, as shown in Equation (37). 

We have completed a large (> 100) galaxy cluster survey in Brownstein and Moffat (2006a) that provides a statis- 
tically significant answer to the question of how much dark matter is expected. The relative abundance across the 
scales of the X-ray clusters would imply that the K-map should peak at ~ 1.0 as opposed to the November 15, 2006 
data release (Clowe et al, 2006b), which peaks at a value of k ~ 0.38 for the main cluster K-map peak. 

Where is the missing dark matter? 

The dark matter paradigm cannot statistically account for such an observation without resorting to further ad-hoc 
explanation, so the question of "missing matter" may be irrelevant. The problem of "extra gravity" due to MOG may 
be a step in the right direction, with our solution as the subject of the next section. 

B. The MOG Solution 

The MOG K-model we have developed in Equations (41) and (42) can now be applied to the Bullet Cluster 1E0657- 
558. In order to model the K-map, shown in Figures 12 and 14a, in MOG, we must integrate the product of the 
dimensionless gravitational coupling, Q{r) of Equation (18) with the mass density, p(r), over the volume of the Bullet 
Cluster 1E0657-558. As discussed in Section III.C, we may integrate the surface density S-map data according to 
Equation (54) to obtain the integrated cluster gas mass about concentric cylinders centered about R = 0. However, 
what we require in the calculation of the Q{r) of Equation (18) is the integrated mass about concentric spheres, 

M(r)= del sin0d</> / r' 2 dr' p(r' ,6,<f>). (71) 
Jo Jo Jo 

The calculation of M(r) using Equation (71) is non-trivial, and we proceed by making a suitable approximation to the 
density, p(r). As discussed in Section II. A, we may apply a sequence of approximations to develop a MOG solution. 

th -order approximation: We have obtained the spherically integrated best-fit King /3-model for the main cluster gas 
mass of Equations (21) and (22) using the parameters listed in Table 2. By neglecting the subcluster, we 
have generated a base-line solution similar to every other spherically symmetric galaxy cluster. The result of 
substituting the spherically symmetric King j3- model mass profile of Equation (22) into Equation (18) is shown 
in Figure 7 for the dimensionless gravitational coupling, Q{r). We obtain a zeroth order, spherically symmetric 
approximation to the K-map by substituting p(r) and Q{r) into Equations (41) and (42) and integrating over 
the volume. 

MOG Center: Treat the subcluster as a perturbation, and shift the origin of Q(r) toward the subcluster (toward the 
center-of-mass of the system). In this approximation, we continue to use the spherically integrated best- fit 
King /3-modcl for the main cluster gas mass of Equations (21) and (22) using the parameters listed in Table 2, 
but we allow the subcluster to perturb (shift) the origin of Q{r) toward the true center-of-mass of the system. 
The integrals of Equations (41) and (42) become nontrivial as the integrand is no longer spherically symmetric. 
We were able to obtain a full numerical integration, but the computation proved to be too time consuming 
(~ 70,000 numerical integrations to cover the 185 x 185 pixel 2 ) to treat by means of a nonlinear least-squares 
fitting routine. 

Projective approximation: Approximate M(r) with the concentric cylinder mass, M{R) of Equation (54), calculated 
directly from the S-map, where the pixel-by-pixel sum proceeds from the MOG center. This is a poor approxi- 
mation for small R (few pixels), but becomes very good for large R (many pixels), as can be seen by comparing 
Figures 8 and 9. The cylindrical mass, M(R), can be computed directly from the E-map data using Equation 
(54) without the need of a King /3-model. 

Isothermal /3-model approximation: Treat the subcluster as a perturbation, and utilize the isothermal /3-model of Equa- 
tion (68) to approximate M(r) and shift that toward the subcluster (toward the MOG center). This is a fully 
analytic expression, allowing ease of integration and an iterative fitting routine. 
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We will base our analysis on our best-fit determined by the isothermal j3- model approximation of Equation (68), 
with G(r) located at the MOG center, a distance away from the main cluster S-map peak toward the subcluster S-map 
peak. Since the dimensionless gravitational coupling, Q{r), of Equation (18) depends on M(r), and the isothermal 
/3-model of Equation (68) depends on Q{r), we must solve the system simultaneously. Let us rewrite Equation (18) 



where 



We may solve Equation (72) for 



(G(r) 1) 

and equate this to the isothermal /3-model of Equation (68) : 

gr) _ M N (r) 
(G(r)-lf G(r) ' 

and we may solve the quadratic equation for 



;i,-)= M^i .-..xi»(-^J (i + L-)\ 2 . (73) 



M(r) = ^ '— , (74) 



(75) 




G(r) = 1 + 7^H^< 1 + \HTT 1 + 1 > ( 76 ) 



where Mn is the isothermal /3-model Newtonian dynamic mass of Equation (65). The result of Equation (76) is a fully 
analytic function, as opposed to a hypergeometric integral, and may be easily computed across the full re-map. There 
was a noticeable parameter degeneracy in choosing the MOG center, the best-fit corresponded to a distance of 140 
kpc away from the main cluster S-map peak toward the subcluster S-map peak. This is reasonable since the S-map 
and re-map data are two-dimensional "surface projections" , due to the line-of-sight integral. The full simulation was 
run iteratively over a range of positions for the MOG center, while covarying the MOG parameters, Mq and ro, the 
MOG mass scale and MOG range parameter, respectively. This yielded a best-fit MOG model for the dimensionless 
gravitational coupling, Q{r) of Equation (18), where r is the distance from the MOG center, as listed in Table 3. Our 
best-fit to the re-map corresponds to the MOG re-model of Equations (41) and (42). The best-fit location of the MOG 
center is provided on the S-map in Figure 3, and the coordinates are listed in Table 1. 

We show a 3D visualization of the convergence re-map data in Figure 14a. The th -order approximation result, 
rescaled is shown in Figure 14b. The twin humps at the peaks of our prediction will be a generic prediction for any 
spherically symmetric galaxy cluster (non- mergers) . The best- fit MOG re-model is shown in Figure 14c along the line 
connecting the MOG center to the main cluster re-map peak. We show a 3D visualization of the full convergence 
re-map model in Figure 14d. 



C. Including the Galaxies 

In considering the MOG re-model resulting from the th -order approximation, as shown in Figure 14b, it is tempting 
to try to explain the entire convergence re-map by the X-ray gas mass, just by shifting the solid-black line to the left 
by ~ 200 kpc, but then there would be no way to explain the subcluster re-map peak. We next proceed to account 
for the effect of the subcluster on the dimensionless gravitational coupling, G(r), of MOG, as shown in the best-fit 
re-model of Figures 14c and 14d. Remarkably, as the MOG center is separated from the main cluster S-map peak, 
say due to the gravitational effect of the subcluster, the centroid naturally shifts toward the re-map peak, and the 
predicted height of the re-map drops. Let us take the difference between the re-map data and our best-fit re-model, to 
see if there is any "missing mass". We hypothesize that the difference can be explained by including the galaxies, 

«( iE , y ) = Sfog) + Ws,I/) > (77) 
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FIG. 13: The galaxy surface density E-map prediction. 

The prediction of the E-map due to the galaxies as computed by the difference between the K-map and our MOG 
K-model, scaled as surface mass density according to Equation (78). E-map and K-map observed peaks are shown 
for comparison. The central dominant (cD) galaxy of the main cluster, the brightest cluster galaxy (BCG) of the 
subcluster, and the MOG predicted gravitational center are shown. J2000 and map (x,y) coordinates are listed 
in Table 1. Component masses (integrated within a 100 kpc radius aperture) for the main and subcluster, the 
MOG center and the total predicted baryonic mass, Mb ary , for the Bullet Cluster 1E0657-558 are shown in Table 5. 



where S is the weighted surface mass density of Equation (42), and the best-fit K-model of S/S c is derived from 
Equations (41) and (42). Therefore, the galaxies contribute a "measurable" surface mass density, 

S ga iax(a;,y) « ^7 \ , (78) 

y{x,y) 

where Q{x, y) corresponds to the best-fit model of Equation (18) listed in Table 3. The result of the galaxy subtraction 
of Equation (78) is shown in Figure 13. Now we may interpret Figure 14d as the total convergence K-map where the 
black surface is the contribution from the weighted surface density of the ICM gas, E/E c , and the red surface is the 
remainder of the K-map due to the contribution of the weighted surface density of the galaxies, S ga i ax /S c . We may 
calculate the total mass of the galaxies, 

Af galax = J T, galax (x' ,y')dx' dy' . (79) 
We were able to perform the integration within a 100 kpc radius aperture about the main cluster cD and subcluster 
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FIG. 14: The convergence K-map November 15, 2006 

The best-fit MOG /-c-model is shown in solid black in 
November 15, 2006 data release (Clowe et at, 2006b) 
14c and 14d. The scaled E-map, E(x,2/)/E c data is 



rfi i — r — r 




(d) 3D visualization of the best-fit K-model 

data release (Clowe er a/., 2006b) and our K-models. 

Figures 14b, 14c and 14d. The convergence K-map 
is shown as Figure 14a and in red in Figures 14b, 
shown in short-dashed green, also shown in Figure 4. 



BCG, separately, the results of which are listed in Table 5, where they are compared with the upper limits on 
galaxy masses set by HST observations. If the hypothesis that the predicted M ga i ax is below the bound set by HST 
observations is true, then it follows that 

Mbary = Mgas + Mgaiax, (80) 

requires no addition of non-baryonic dark matter. The results of our best-fit for Af gas , M ga i ax and Mb ary of Equation 
(80) are listed in Table 5. 
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Component 


Main cluster 


Subcluster 


Central ICM 


Total 


M gas 

Mgalax 


7.0 x 10 12 Mq 
1.8 x 10 12 Mq 


5.8 x 10 12 Mq 
3.1 x 10 12 Mq 


6.3 x 10 12 Mq 

2.4 x 10 10 Mq 


2.2 x 10 14 M Q 
3.8 x 10 13 M© 


-^bary 


8.8 x 10 12 M 


9.0 x 10 12 Mq 


4.9 x 10 12 Mq 


2.6 x 10 14 M 


Mum 


2.1 x 10 13 Mq 


1.7 x 10 13 Mq 


1.4 x 10 13 Mq 


6.8 x 10 14 M Q 


M ga i ax / Mg m 
M sas /Mum 


26% 
33% 


53% 
34% 


0.4% 
45% 


17% 
32% 



TABLE 5: Summary 
of component mass predictions. 

Component masses (integrated 
within a 100 kpc radius aper- 
ture) for the main and subclus- 
ter and the MOG center. The 
total predicted mass for the Bul- 
let Cluster 1E0657-558 is in- 
tegrated over the full E-map. 



D. Dark Matter 

From the alternative point-of-view, dark matter is hypothesized to account for all of the "missing mass" which 
results in applying Newton/Einstein gravity. This means, for the November 15, 2006 data release (Bradac et al., 
2006; Clowe et al., 2006a,b,c), that the "detected" dark matter must contribute a surface mass density, 

£dm(e,#) ~ K(x,y)T, c - T,(x,y), (81) 

with an associated total mass, 



Mum =1 Xum(x', y')dx'dy'. (82) 
The integral of Equation (82) becomes trivial upon substitution of Equation (81): 

Mum = S c / k(x', y')dx' dy' - [ Z(x',y')dx' dy' (83) 



where we have neglected M ga i ax in Equation (83), because as is usually argued, the contribution from the galaxies 
in the dark matter paradigm is < 1 - 4% of M to tai- The calculation of Mum in Equation (83) was performed by a 
pixel-by-pixel sum over the convergence K-map data and surface density S-map data, within a 100 kpc radius aperture 
around the main and subcluster K-map peaks, respectively. The result of our computation is included in Table 5. We 
emphasize, here, that for each of the main cluster, subcluster and total S-map, our results of Table 5 indicate that 

Mbary = Mg as + M ga l ax < M DM (84) 

implying that we have successfully put Bullet Cluster 1E0657-558 on a lean diet! This seems to us to be a proper 
use of Occam's razor. The mass ratios, Af ga i ax /M gas , for the main and subcluster and central ICM are shown at the 
bottom of Table 5. The result of M ga i ax /M gas w 0.4% in the central ICM is due to the excellent fit in MOG across the 
hundreds of kpc separating the main and subcluster. The dark matter result of M gas /MuM ~ 45% in the central ICM 
implies that the evolutionary scenario does not lead to a spatial dissociation between the dark matter and the ICM gas, 
which indicates that the merger is ongoing. In contrast, the MOG result shows a true dissociation between the galaxies 
and the ICM gas as required by the evolutionary scenario. The baryon to dark matter fraction over the full S-map 
is 32%, which is significantly higher than the A-CDM cosmological baryon mass-fraction of 171^2% (Spergel et al, 
2007). The distribution of mass predicted by MOG vs. the dark matter paradigm is shown in Figure 15. 



V. CONCLUSIONS 

The modified gravity (MOG) theory provides a fit to the K-map of the November 15, 2006 data release (Clowe et at, 
2006b). The model, derived purely from the X-ray imaging S-map observations combined with the galaxy S-map 
predicted by MOG, accounts for the K-map peak amplitudes and their spatial dissociation without the introduction 
of nonbaryonic dark matter. 

The question of the internal degrees of freedom in MOG has to be further investigated. It would be desirable to 
derive a theoretical prediction from the MOG field equations that fits the empirically determined mass and distance 
scales in Figure 2. 

It could be argued that any modified theory designed to solve the dark matter conundrum, such as MOND or 
MOG, has less freedom than dark matter. So the important question to resolve is precisely how much freedom the 
MOG solution has. On one hand we said, definitely, that there was no freedom in choosing the pair of Mq and ro 
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(a) Scaled surface density for the MOG (b) Scaled surface 

predicted galaxies, E ga iax/E c , and the MOG predicted density of dark matter, Edm/E c , compared to ICM gas. 

visible baryons, Eb a ry/E c , compared to the ICM gas. 

FIG. 15: Plot of the scaled surface density E/E c along the line connecting the main cluster E-map peak with the main cD. 

In Figure 15a, the prediction of Equation (78) for the galaxies is shown in long-dashed magenta, and the prediction of Equation 
(80) for the visible baryonic mass is shown in solid brown. The calculation of Equation (81) for dark matter is shown in Figure 
15b in dash-dot black. The ICM gas distribution inferred from the E-map data is shown in short-dashed green on each plot. 

for the main cluster since it was well described by an isothermal sphere, to an excellent approximation. We further 
argued that the subcluster was (per mass) a small perturbation to the ICM. But if MOG has more freedom than 
MOND, but less freedom than Dark Matter, then what is the additional degree of freedom that enters the Bullet 
Cluster observations? 

The question is resolved in that there is a physical degree of freedom due to a lack of spherical symmetry in the 
Bullet Cluster, and whence the galaxies sped outward, beyond the ICM gas clouds which lagged behind - effectively 
allowing the galaxies to climb out of the spherical minimum of the Newtonian core where MOG effects are small 
(inside the MOG range ro) upwards along the divergence of the stress-energy tensor (Newtonian potential, if you 
prefer a simple choice) towards the far infrared region of large gravitational coupling, Goo. 

In fact, the Bullet Cluster data results describe, to a remarkable precision, a simple King /3-Model. Our analysis, 
with the result to the best-fit shown in Table 2, uniquely determines the mass profile p(r) of Equation (21) used 
throughout our computations. We permitted only a single further degree of freedom to account for the fits of Figure 
14c and the predictions of Figures 13 and 15; this was the location of the MOG center, where the gravitational 
coupling, (?(0) — > 1, is a minimum at the Newtonian core. Remarkably, the data did not permit a vanishing MOG 
center, with respect to the peak of the ICM gas p(0). We have shown the location of the MOG center as determined 
by a numerical simulation of convergence map according to Equations (41) and (42) in each of Figures 3, 10, 12 and 
13 and provided the coordinates in Table 1. 

The surface density E-map derived from X-ray imaging observations is separable into the main cluster and the 
subcluster subtracted surface density S-map through a low % 2 -fitting King /3-model. Following the (> 100) galaxy 
cluster survey of Brownstein and Moffat (2006a), we have derived a parameter-free (unique) prediction for the X-ray 
temperature of the Bullet Cluster 1E0657-558 which has already been experimentally confirmed. In Equations (41) 
and (42), we have derived a weighted surface mass density, S, from the convergence K-map which produced a best-fit 
model (Figures 14c and 14d and Table 3). We have computed the dark matter and the MOG predicted galaxies and 
baryons (Figure 15), and noted the tremendous predictive power of MOG as a means of utilizing strong and weak 
gravitational lensing to do galactic photometry - a powerful tool simply not provided by any candidate dark matter 
(Figure 13). The predictions for galaxy photometry will be the subject of future investigations in MOG, and the 
availability of weak and strong gravitational lensing surveys will prove invaluable in the future. 

Although dark matter allows us to continue to use Einstein (weak-field) and Newtonian gravity theory, these theories 
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may be misleading at astro-physical scales. By searching for dark matter, we may have arrived at a means to answer 
one of the most fundamental questions remaining in astrophysics and cosmology: How much matter (energy) is there 
in the Universe and how is it distributed? For the Bullet Cluster 1E0657-558, dark matter dominance is a ready 
answer, but in MOG we must answer the question with only the visible contribution of galaxies, ICM gas and gravity. 
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